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ABSTRACT 

The present observational understanding of the evolution of the mass accretion rates (Mace) in pre- 
main sequence stars is limited by the lack of accurate measurements of Mace over homogeneous and 
large statistical samples of young stars. Such observational effort is needed to properly constrain the 
theory of star formation and disk evolution. Based on HST/WFPC2 observations, we present a study 
of Mace for a sample of ^ 700 sources in the Orion Nebula Cluster, ranging from the Hydrogen-burning 
limit to M^, ~ 2Mq. We derive Mace from both the U-hand excess and the Ha luminosity {Lna), after 
determining empirically both the shape of the typical accretion spectrum across the Balmer jump and 
the relation between the accretion luminosity (iaec) and Lh^i that is L^cc/Lq = (1.31±O.O3)-L/fc[/^0 
+ (2.63± 0.13). Given our large statistical sample, we are able to accurately investigate relations 
between Mace and the parameters of the central star such as mass and age. We clearly find Mace to 
increase with stellar mass, and decrease over evolutionary time, but we also find strong evidence that 
the decay of Mace with stellar age occurs over longer timescales for more massive PMS stars. Our best 
fit relation between these parameters is given by: log(Maec/M0-yr)=(-5.12 ± 0.86) -(0.46 ± 0.13) 
•log(Vyr) -(5.75 ± 1.47)-log(M*/MQ) + (1.17 ± 0.23)- log(t/yr) • log(M*/M0). These results also 
suggest that the similarity solution model could be revised for sources with > 0.5Mq. Finally, we 
do not find a clear trend indicating environmental effects on the accretion properties of the sources. 



1. INTRODUCTION 

During the Pre-Main-Sequence (PMS) phase of stel- 
lar evolution, the interaction between the forming star 
and the surrounding disk is regulated by the accretion of 
disk m aterial along the field li nes of the stellar magneto- 



sphere ( Hartmann et aL|l998 ). The gravitational energy 
released as the material falls along accretion columns and 
hits th e stellar surface cr eates a characteristic shock spec- 
trum (Calvet & Gullbring 1998^), with excess emission 
especially str ong in the~Balmer con t inuum and recombi- 
natio n lines ( Gullbring et al.|[l998 [ Calvet & Gullbring 
19981. The relative accretioii luminosity (j^aec ) can be 



measured with spectros copic (Valenti et al.1'1993 Her 



czeg & Hi llenbrand|2008 ) o r photometric ( Gullbring et al. 
Tgg8| iRobberto et al.|2TO4l|De Marchi et al.|2010| ) meth- 
ods. The mass accre tion rate, Maeci is the n estimated 
through the relation ( Hartmann et al.||1998 ): 



Ma 



0.8G'M, ' 



(1) 



where i?* and M* are the radius and the mass of the 
star, respectively, and the factor 0.8 accounts for the as- 
sumption that the infall origin ates at a magnetospheric 
radius i?^ =5 ( |Shu et al.||1994[ ). 
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The mass accretion rate generally decreases over time 
during the first few Myr of PMS stellar evolution, as the 
circumstellar disks disp erse their gaseous content on a 
timescale of ^ 3-5 Myr ( [Haisch et aT]|2001[ |Dahm||2005" 
Fedele et al.||2010 ). Moreover, Mace is expected to scale 
with stellar mass. The evolution of Mace vs. time as a 
function of the stellar and disk mass represen ts a key as- 
pect o f PMS evolution and planet formation. |Hartmanii| 



et al. 



( 1998 ) showed that while typical Mace values, at 



any given stellar age, can differ star by star of up to about 
two orders of magnitudes, on average they decrease ex- 
ponentially with the stellar age t, i.e. MaecOc i"'' with 
77 ~ 1.5. Several studies, targeting different young stellar 
clusters, hav e conf irmed this general trend (e.g. Rob^ 
berto et aL||2004[ |Sicilia-Aguilar et al. 2010), but the 



uncertainties m the age es timate ( Hartmann et al.^^1997 
Baraffe & Chabrierl 2010) and the scarcity ot ricn anc 



homogeneous samples limit the accurate assessment of 
this dependence. 
For what concerns th e stellar mass dependence of M., 

Mf with b 



Muzerolle et al. 



(2003) found that Ma 



acc; 

2, 



although recent studies (e.g., Rigliaco et al. 2011 1 suggest 
smaller values of b. Due to the large spread in the Aface 
values for a given (up to two orders of magnitudes, 
Rigliaco et al.||20TT ), this second relation has also never 



been accurately constrained. Although observational un- 
certainties and intrinsic variability will always contribute 
in scattering the measured mass accretion rates, accurate 
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measurements of Mace on a large sample of PMS stars 
may allow us to assess more precisely the value of the two 
power law exponents (77 and b). This is a critical step, 
as they can be tied to the theoretical model for the disk 
evolution and structure, providing unique constraints to 
the initial conditio ns of pla net formation (e.g. , Alexander, 
fc Armitage]|2006[ lLodato||2008[ ) . 



As the nearest (d ~ 414 ± 7 pc, Menten et al. 20071 
site of massive star formation, the Orion INebula Cluster 
(ONC) provides a standard benchmark for star forma- 
tion studies. The ONC population has been studied in 
depth, and determination of the individual stellar pa- 
rameters are available for a larg e fraction of members 



In particular |H illcnbrand| ( |1997[ ), |Da Rio et aT] ( |2010 l 
and Da Rio et al.| (j201^have used both spectroscopic 



and photometric techniques to derive the spectral type 
of more than 1700 ONC sources. The corresponding ini- 
tial mass function, ranging from the Brown Dwarfs ( BD) 
regime to the 6 star 9^0n-C, peaks at 0.2-0.3 Mq (Da 
Rio et al.|2012 ), while the cluster mean age is ~ 2.2 Myr 
with evidence f or an a ge spread of the order of ^ 2 Myr 



( Reggiani et al.pOll ). Given the wealth of information 
on the individual cluster members, the ONC is ideally 
suited to conduct an extensive study of the mass accre- 
tion process in PMS stars. 

In this paper we present the results of such a study 
based on the HST/WFPC2 survey of the Orion Nebula 
Cluster (GO 10246, P.I. M. Robberto). Both U-band and 

Ha data are used to estimate Mace for ^ 700 sources 
and to analyze the relations between this parameter and 
the main stellar parameters (age, M*). In Sec. [2] we 
illustrate the observation, data reduction and analysis 
of our WFPC2 photometry. Sec. [3] presents our derived 
modeling of the photospheric colors for our sources and 
for the typical accretion spectrum, and the methods used 
to obtain the Lace based on our data, both using our UBI 
Diagram method and from the Ha photometry, while in 
UJwe present the derived quantities (Lace, Mace) obtained 
Doth from the measurement of the J7-band or from Lhu 
in Sec. [5] we study the evolution of these parameters as 
a function the main stellar parameter. Finally, in Sec. |6] 
we summarize our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observations 

In this work we consider the Hubble Space Telescope 
(HST) observations obtained for the HST Treasury Pro- 
gram on the ONC ( [Robberto et al.||2012[ ). In particular, 
we use the WFPC2 observations carried out between Oc- 
tober 2004 and April 2005 in the filters F336W, F439W, 
F656N and F814W. One of the advantages of this data 
set is the nealry simultaneous (within less than one hour) 
imaging obtained in subsequent exposures in U (F336W) 
and Ha (F656N) bands, which, as will be described in 
the next sections, both allow an estimate of the accretion 
l uminosity. 

Robberto et al. ( 2012) details the data acquisition and 



reduction strategy of all the photometry obtained within 
the Treasury Program; here we summarize only the key 
aspects of the WFPC2 dataset, which are most relevant 
to our study. The F336W filter is roughly analogous to 
the standard Stromgren u filter, with a long wavelength 
cutoff at A ~ 3600 A, shortward the Balmer jump. This 



WFPC2 filter, however, is affected by "red leak" , which 
is a residual transmission at longer wavelengths (A ~ 
7300 A). This effect alters the measured [/—band fluxes 
for red sources, and in Appendix [A| we explain how the 
correction for this leak is performed. We hereafter refer 
to the corrected U band as U, using the notation Unoieak 
only when needed, for clarity. The F439W filter closely 
matches the standard Johnson B passband, while the 
F656N filter is a narrow-band filter centered on the Ha 
line, narrow enough to exclude the adjacent [Nil] doublet 
lines. The last filter, F814W, is similar to the Cousin Ic 
filter; it is fundamental for our study both to determine 
the red leak contamination of the [/-band, and to derive 
the bolometric luminosity of the ONC members. The 
exposure times for each field are 800 s in F336W, 80 s in 
F439W, 400 s in F656N and 10 s exposure in F814W. 

Aperture photometry was obtained by choosing aper- 
ture radii of O'.'l and (X'S, corresponding to 5 and 11 pixels 
on the WFl (PC) and 2 and 5 pixels on WF2-4 chips. 
The sky annulus has always been taken between 1" and 
I'.'S, corresponding to 20 and 30 pixels with WFl and 
10 and 15 pixels with WF2-4. At the time the obser- 
vations were carried out, the WFPC2 instrument had 
already accumulated about 11 years of total radiation 
dose in space environment, and therefore was affected by 
non-negligible Charge Transfer Efficiency (GTE) losses. 
The brightness of the Orion Nebula background, espe- 
cially in our four broad-band filters and at the center of 
the region, mitigates the problem but still we had to ap- 
ply a CTE correction to th e measured counts following 
the recipe of Dolphin (2009). We estimated CTE correc- 
tion errors using a Monte Carlo propagation of the errors 
terms in the equation of .Dolphin (2009). 

The final photometric errors have~been derived by 
adding in quadrature the errors associated to the mea- 
sured counts, zero point, CTE correction, aperture cor- 
rection to the infinity and, in the case of the O'.'l aper- 
tures, an extra aperture correction to the 0'.'5 radius. 
This last term is usually dominating, ensuring that the 
0'.'5 photometry is generally more accurate and therefore 
preferred, except for the faintest sources detected in very 
few pixels. 

2.2. Source selection 

Our WFPC2 photometric catalog includes 1643 
sources in total with detection in at least one of the 
four bands. Among t hese, we have rem oved 27 known 
proplyds, identified by Ricci et al. (2008) from the ACS 
imaging, since the fluxes of these are potentially contami- 
nated by circumstellar emission. We have also eliminated 
144 close binary systems (d < 3 pix ^ 75-100 AU), re- 
solved either in the ACS or the WFPC2 imaging, where 
the spectral type classification is more uncertain. For 
well separated visual binaries, we carefully inspected the 
individual sources to rule out possible matching errors 
between the photometry in different bands. 

In the remaining sample, C/-band photometry is avail- 
able for 1021 stars (^ 60% of the total), among these, 
897 sources (^ 55%) have been detected in all [/, B and 
/ bands. Finally, for 1342 sources (~ 80%) Ha photom- 
etry is available. 
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Figure 1. Our U BI 2CD; the thick hne is our calibrated isochrone 
vahd for the ONC. The arrows indicate two reddening vectors, 
corresponding to Av =2, for two different values of Ry assuming 
the extinction law of |Cardelli et alT| ( |1989[ l 

3. METHOD 



We follow the approach introduced by |Da Rio et al. 
( 2010 1 to obtain simultaneously an estimate of extinction 
(Av) and accretion (Lacc/itot, where Ltot = + ^acc) 
for individual sources from the photometry in 3 optical 
bands together with the T^g of the source. The procedure 
assumes that on a color-color diagram the displacement 
of the observed sources from the theoretical isochrone is 
due to a combination of extinction and accretion. First 
of all accretion excess causes a shift in the colors to- 
wards bluer values, then extinction moves the resulting 
color along the reddening direction. In particular, in a 
two-color diagram, accretion displaces the sources along 
curves that start from the isochrone and converge to a 
point defined by the colors of the pure accretion spec- 
trum. If Toff of a star is known, e.g., from spectroscopy, 
the combination of Ay and ^acc/^tot necessary to repro- 
duce thfi_observBd cofors is unequivocally found. 



In Da Rio et al.| (20101) a BVI two-color diagram was 
used; here instead we construct the J/noicak — B versus 
B — I two-color diagram. Since [/-band is far more sen- 
sitive to accretion than B, this allows us to measure 
L^rr J Ltnt wlth & si guificantly higher accuracy than in 



Da Rio et al. (2010). We will refer from now on to this 
two-color diagram as the 2CD. 

3.1. The Two-Colors Diagram: Extinction and 
Accretion Luminosity 

We con sider the measured T pff for the ONC mem- 
bers from Da Rio et al. (20121, and limit ourselves to 
the sources detected in all the three WFPC2 bands and 
where the error on the Unoieak values is ^0.3 magj^ Our 
sample is then composed by 339 ONC members/~Fig. [T] 
shows the observed 2CD of this sample; the solid line 

^ The phot ometrical errors in all the bands {U,B,I), calculated 
as in Sec. |2.1[ are always lower than 0.2 mag, but the correction for 
the red leak (see Appendix|A]l raises the final error on the measure 
to values up to ~0.29 mag m 15 objects. 



represents our calibrated isochrone (see Sec. 3.1.1 for de- 
tails). 

To model the ex tinction, we assume the reddening law 
1989) with typical galactic reddening 
3.1. Some authors have sug gested a 
larger value of Ry = 5.5 for the ONC (e.g., Costero 



of Cardelli et al 
parameter Ky 



& Peimbert 



Da Rio et a 



1970) and we show this as well in l^'ig. ^ 
^2010), however, found Ry = 3.1 to be 



more appropriate, on average, for the ONC population. 
Therefore, we will adopt this value for our analysis. 

3.1.1. Color Calibration 

As required by our method, we need to assume an 
isochrone that defines the intrinsic, photospheric, colors 
as a function of Tp g , for th e ONC sources in our 2CD. 
As shown by Da Rio et al. (2010), current atmosphere 
models are still unable to lead to accurate predictions of 
optical colors for young (~ 2 Myr) cool stars, and em- 
pirical calibrations of photospheric colors have not been 
specifically carried out for young PMS stars. Also, in 
that work the authors show that, for the same Tes, the 
photospheric colors of ^ 2 Myr PMS sources are not in 
agreement with those of MS dwarfs. Thus we cannot ei- 
ther rely on calibrations of colors valid for more evolved 
populations. This problem therefore requires a particular 
care. We explain in Appendix IB] the procedure used to 
obtain the correct isochrone calioration for our 2CD, and 
we report in Table [T] the colors of the isochrone. It could 
be argued that the assumption of a single isochrone in 
our 2CD, given a possible age spread in Orion, may lead 
to false predictions. This is not the case, for 2 reasons: 
a) the 2CD is, by its nature, independent on stellar lu- 
minosities or radii; therefore age variations for the same 
Toff can result in differences in the photospheric colors 
only due to changes in the stellar surface gravity; b) in 
the optical, the difference in colors between PMS and MS 
stars ca n be of up to several tenths of a magnitude (Da 
Rio et al. 2010). This offset , scaled down to the modest 
age spread within the ONC ([Reggiani et al. 2011 ), results 
in negligible differences in the colors vs. les relations. 

In Fig.fll we present our calibrated isochrone, together 
with the observed sources. If ongoing mass accretion did 
not alter our observed colors, the sources would have 
been located on the model line, with only a spread given 
by extinction. Instead, a number of ONC members are 
displaced towards bluer colors, as expected by the pres- 
ence of accretion flux excess. 

For our subsequent analysis, and specifically for the 
derivation of bolometric stellar luminosities, we also re- 
quire the bolometric corrections {BC) in the /—band as 
a function of Tcff. These are defined as the difference 
between the bolometric magnitude of a star and the ap- 
parent /—band magnitude. 



: 2.5 log 



BC = Mboi - Mx 



IhjT) . , 

. Jh{T)-dX I Jh{Q)-dX 



(2) 



We computed the BCs by perform ing synthetic pho- 
tometry on the spectra of Allard et al. (2010), and cor- 
rected the results based oii the same offsets we determi- 
ned in our calibration of the colors, to correct the Tband 
intrinsic magnitudes. The values of this empirical cal- 
ibrated bolometric correction for various spectral types 
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Table 1 

Intrinsic colors for the WFPC2 bands F336W 
(U), F439W (B) and F814W (I) and Bolometric 
Correction f or the F814W (I) ba nd obtained 
assumins; [Allard et aT] ||2010[| models. 



Spectral 


Tcff 

(K) 


{U-B) 


(B-I) 


(BCi) 


Type 


(mag) 


(mag) 


(mag) 


A5 


8200 


0.126 


0.285 


-0.450 


A7 


7850 


0.106 


0.424 


-0.365 


FO 


7200 


0.066 


0.673 


-0.236 


F2 


6890 


0.017 


0.809 


-0.180 


F5 


6440 


-0.013 


1.022 


-0.105 


F8 


6200 


-0.002 


1.151 


-0.068 


GO 


6030 


0.022 


1.247 


-0.044 


G2 


5860 


0.054 


1.352 


-0.021 


G5 


5770 


0.075 


1.406 


-0.012 


G8 


5570 


0.142 


1.530 


0.007 


Kl 


5080 


0.457 


1.871 


0.032 


K2 


4900 


0.597 


2.024 


0.033 


K3 


4730 


0.710 


2.188 


0.027 


K4 


4590 


0.790 


2.342 


0.017 


K5 


4350 


0.945 


2.636 


-0.010 


K7 


4060 


1.111 


2.997 


-0.055 


MO 


3850 


1.234 


3.272 


-0.108 


Ml 


3705 


1.313 


3.460 


-0.159 


M2 


3560 


1.412 


3.646 


-0.218 


M3 


3415 


1.543 


3.890 


-0.288 


M4 


3270 


1.665 


4.266 


-0.378 


M5 


3125 


1.771 


4.722 


-0.515 


M6 


2990 


1.771 


5.143 


-0.762 


M7 


2880 


1.734 


5.513 


-1.050 



are also shown in Table [T] The BC/ varies between ~ 
0.03 and -0.70 magnitudes in the range of temperatures 
of our sources (3000 K < Teff < 8000 K). 
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Figure 2. 2CD of our ONC population. Each source has been cor- 
rected for extinction and represented as a normalized 2D gaussian, 
corresponding to the photometric errors. The sources arc color- 
coded according to their T^fi, as shown in the scale on the right 
of the plot. The thick line represents our calibrated isochrone, for 
no accretion. The thin lines represent the simulated displacements 
(for different T^g) from the photospheric colors obtained by adding 
an increasing amount of Lacc/^tot, assuming the accretion spec- 
trum analogous to that of CG98. Dashed lines indicate constant 
ratios Li^cc/Ltot 



3.1.2. Accretion spectrum 
Spectra of the f lux excess produced by acc retion have 



been modeled by Calvet & Gullbring (11998) (hereafter 
CG98) using a hydrodynamic treatment for the infalling 
material and heated photosphere. It is possible, however, 
to reproduce a typical accretion spectrum using sim- 
pler recipes, for instance modelling the radiati ve transfer 
within a slab of dense gas ( Valenti et al.|1993 1 or assum- 
ing suitable combination ot an optically thm emission 
generated in the preshock region and an optically thick 
emission generated by the heated photosphere. In the 
latter case the two components can be added with a rel- 
ative fraction of a bout 1/4 and 3/4, respectively ( |Gull- 
bring et al. 2000). For simplicity, we follow this latter 
approach; tor the optically thick component, we consider 
a 8000 K black body, and for the optically thin one, we 
modeled the emerging spectrum from a slab of gas using 
the version 10.00 of the Cl oudy software, a pho toioniza- 
tion code last described by Ferland et al. (19981. Specif- 



10^ 



ically, we assumed a slab with density 
and temperature of ~ 20000 K. The particular choice of 
these two parameters has been selected after many tri- 
als, in order to obtain an accretion spectrum as similar 
as possible to that of CG98. For the far-UV par t of the 
spectrum (A < 3000 A) we follow the results of France 
et al. (2011), and assume a linear decrease of hux at 
wavelengths shortward of the Balmer continuum. 

We use our photometrical colors and produce a red- 
dening co rrected 2CD (F i g. [2| assuming again the Ay 
values of Da Rio et al. (2012). In this way, the dis- 



placements ot the UiN C members from the isochrone are 



solely due to accretion and we can investigate this as- 
pect, neglecting the extinction contribution. In this plot 
sources are represented as equally normalized 2D gaus- 
sian, with widths corresponding to the individual pho- 
tometric error. This facilitates our visual inspection, 
because, in this representation, sources with small er- 
rors appear sharp and bright while sources with poorer 
photometry are broader and less luminous, contribut- 
ing to the plot mostly through their relative number. 
The sources are also color-coded according to their Tgff, 
highlighting that the dereddened sources are generally 
located at bluer colors than the photospheric ones, along 
stripes that depend on Tcff. 

In Fig. [2] we also show the predicted displacement from 
the photospheric isochrone, obtained by adding a com- 
ponent of the basic accretion spectrum to the theoretical 
(cahbrated) spectra. Specifically, increasing iacc/^tot 
from to 1, the colors are moved along lines that start 
from the theoretical isochrone, at any stellar Tcff, and 
converge to the colors of the pure accretion spectrum, lo- 
cated on lower left end of the plot. For this computation, 
we have used the basic accretion spectrum, analogous to 
that of CG98. It is evident that, according to this model, 
accretion affects both the {B — I) and the {U — B) col- 
ors, and this prediction is somewhat compatible with the 
extinction-corrected photometry of some ONC members. 
On the other hand, there are many sources observed at 
{U — B) < —1, a region of this diagram which is not 
reachable according to this accretion spectrum model. 
This suggests that for these stars the Balmer jump must 
be higher than that of our accretion spectrum. 
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Figure 3. Same as Fig. [2] but adding to the stellar spectra a sum 
of the accretion spectrum produced used for the model of Fig. [2] 
and a pure HII region with decreasing ratio of this component to 
the first for increasing Lace. 

To correct for the discrepancy we have included an- 
other component to this spectrum in order to have a 
better agreement with our data. This component is mod- 
elled with an HII region spectrum, also simulated with 
Cloudy, with the same density as the one previously 
assumed for the optically thin part but a much lower 
temperature (~ 3000 K). This is added to the accretion 
spectrum in such a way that it becomes relevant only for 
small amounts of mass accretion. The result is shown 
m Fig.^ the new model predicts, consistently with our 
data, a larger displacement in {U — B) for low values of 
■t'acc/-t'tot- For our subsequent analysis we will therefore 
consider this approach to model the color displacement 
due to mass accretion. 

3.1.3. Final Two-Colors Diagram 

From the modeling shown in Fig.pl we derive iacc/^tot 
for all our ONC sou rces. For eacn s ource , we consider 
its known Teg from Da Rio et al. (2012), and there- 



fore the unique accretion displacement curve in the 2CD 
that starts from isochrone point corresponding to the 
stellar temperature of the source. We then consider 
the observed photometry of the source in the 2CD, and 
deredden it along the reddening vector until a intercep- 
tion with the accretion curve is found. This intersec- 
tion defines, geometrically, both the amount of Ay and 
iacc/itot for the star. 

Using this method, we obtain a solution for 271 
sources. In order to quantify the uncertainty in the re- 
sults, we use a Monte Carlo method: for each source we 
displace randomly the photometry in each band within 
the photometric errors, assumed gaussian. We also ran- 
domly change the stellar Tgff accounting for the uncer- 
tainty in the spectral types from spect roscopy, equal to 
±1 spectral subtype ( Hillenbrand||1997 ) . We iterate this 
approach 1000 times for each star, deriving each time Ay 



and iacc/itot values. From the 1000 results of Ay and 
-^acc/itot for each source, we determine the 1, 2 and 3 a 
intervals of the distribution of results. If an intersection 
is found for all the 1000 iterations, we assume our result 
to have a 3 sigma confidence; on the other hand if for 
some of the iterations an intersection is not found (be- 
cause, e.g., the colors lie outside the color range reachable 
by the accretion model plus extinction), a lower confi- 
dence level is associated to the result. We report these 
values in Table [2] Neglecting sources with confidence 
lower than 1 cr, our sample comprises 245 sources. In 
Fig. |4]we show the distribution of Ay and Lacc/Ltot for 
this sample. Our extinctions vary from ~0 to ~ 5, with 
a mean value of ^ 0.98 ± 0.06 magnitudes and only 7 
sources show a negative value of Ay, the lowest being 
Ay — — 0.18± 0.03 mag. We assign to these sources Ay 
= and retain them in our a nalysis. Comp aring our 



results with those obtained by Da Rio et al. ^2012 ) we 
obtain a very good match, measuring of a distribution of 
differences AAy = Ay^^.^ ^^^^ — ^vbaRio peaked at -0.03 
mag with a sigma of 0.31 mag. This testifies the accuracy 
of our method. 

From now on we will use these values of Ay and Lace 
for all the sources for which this method gave a result. On 
the other hand, a small difference in the estimate of Lace 



with Da Rio et al. (2010) is expected, given the higher 
sensitivity of our 2Cd to this quantity with respect to 
their analysis. 

The Monte Carlo approach allows to estimate the er- 
rors on our parameters. We derive individual errors from 
the standard deviation of the values obtained; these will 
be used in our subsequent analysis to estimate the errors 
of the Mace- 

3.2. Ha Luminosity 

Besides the U—hand excess, magnetospheric infall ac- 
cretion processes lead to line emission, in particular the 



hydro gen recombination lines (see e.g., Muzerolle et al 



2001|. Specifically, the Ha excess is often used as an 
indicator of mass accretion, and its intensity can be em- 
pirically related to Lace- 

As anticipated in Sec. |2.2[ our WFPC2 catalog includes 
1342 sources with measured Ha photometry. In order 
to correctly estimate the Equivalent Width of the Ha 
emission (EWhq) one requires an estimate of the pho- 
tospheric continuum at its wavelength. This generally 
depends depends on a nearby broad band flux (e.g., the 
/-band luminosity), as well as both Tgff and Ay of the 
sources. To this pu rpose, we consider temperatures from 



Da Rio et al. (2012) and exti nctions fr om either our 2CD 
analysis (Sec. 3.1.3 1, or from Da Rio et al. (2012) for the 
sources with no results from the 2CD. These required 
parameters reduce our sample to a total of 1026 sources. 
The photospheric continuum is then estimated by means 
for synthetic photometry, computing the photospheric 
(Ha — I) color for each star given its Tes, and subtract- 
ing it from the dcreddened / magnitudes. The Ha ex- 
cess, measured from the difference between the observed 
Ha flux and the expected photospheric flux, is then con- 
verted in terms on equivalent width, or in terms of frac- 
tion of total stellar luminosities. We also stress that the 
WFPC2 F656N filter is broad enough (FWHM-28.5 A~ 
1300 km/s) to contain the entire Ha excess, since this 
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L 
-4 



log L„c./ 



Figure 4. Histograms of the values of Ay and Lacc/Ltot estimated for our sources from the analysis of the 2CD. 



line has a t ypical broadening of ^200-700 km s ^ (^4.5 



13A) (e.g., White fc Basri 



2003 



Mortier et al. 



2011), 



much smaller than the filter width. 

From all the derived values, we consider in our analy- 
sis only sources with 3 < |EW//q,| < 1000. This choice 
allows us to exclude non-accreting sources, or members 
showing a suspiciously high excess, probably due to un- 
certainties or contamination from circumstellar material. 
We also exclude from our analysis sources that are known 
from visual inspection on the HST images to be proplyds, 
binaries or edge-on disks, as already explained in Sec. |2.2[ 
Ultimately, the final sample of sources with available Ha 
luminosity includes 682 ONC members. We report these 
results in Table [3l 

3.3. Stellar Luminosity 

In order to convert the obtained values of Lacc/Ltot j 
an estimate of the stellar luminosity of our sources is re- 
quired. We derive the bolometric luminosity (L* ) of each 
source from the observed / magnitudes, corrected f or the 
effects of extinction and accretion derived in Sec. 13.1.31 
In particular, accretion is taken into account by com- 
puting the excess in magnitudes A/acc in the /—band 
due to the derived amount of accretion for each source. 
These values A/acc (reported in Table [2]) have a mean 
of ~ —0.04 mag, a rather small v alue, since / band i s 
weakly affected by accretion (e.g. Fischer et al. 2011). 
Thus, L» is then derived from: 



log 



0.4 • [M 



= 0.4 • [M,„^p^,,0 



bol,0 



M 



bol, 



"-SC/„ppc2,0 — BCi^ 



(3) 



= 0.4 • [M 



V,0 

^/WFPC2 



{V-L 



WFPC2j0 



/WFPC2 + 

DM], 



where Iwfpc2 refers to magnitudes in the WFPC2 pho- 
tometric system, V indicates the T^-band magnitude in 
the standard J ohnson photometry. We as sume the value 
Mv,Q = 4.83 ( |Binney &: Tre mainc' 2008 ), checking that 
using as a r eference tor the sola r spectrum the synthetic 
spectrum of |Allard et aT] ( |2010p with Teff,0=578O K and 
log 50=4.43 we obtam the same result, and we compute 
through synthetic photometry {V — Iwfpc2)q = 0.70. 
We assume a distan ce for the ONC of d = 414 ± 7 pc 
(Menten et al. 2007), corresponding to a distance mod- 



V 



ulus DM — 8.085. Extinction is conv erted from A 
toAfjyy^^j using the reddening law of Cardelli et al. 
(19891; also here the values of Ay adopted are those 



et al 
FrniH 



obtained wi th the 2CD, if available, and from |Da Rio| 
(2012), when there are no results from the 2CD. 



(Tc ff) are the I-band bolometric correc- 
tions, derived in Sec. |3.1.1[ The derived L* are reported 
in Table [2] in the cases where Ay is estimated with the 
2CD and in Table [3] in the others. 

Errors in are also derived with a Monte Carlo 
procedure, similarly to Sec. |3.1.3| propagating the es- 
timated uncertainties of the photometry, extinction, DM 
and A/acc- 

4. RESULTS: STELLAR PROPERTIES 
4.1. Accretion Luminosity 
4.1.1. iacc from the Two-Colors Diagram 

We conv ert the derived /vacc/^tot from the 2CD 
(Sec. 3.1.3) in terms of Lacc//^0- This is achieved con- 



sidering the stellar luminosities L* derived in Sec. 3.3 
from the relation: 



t 



Lq 1 — {Lacc/Ltot) 



L^ 



(4) 



In Table [4] we report these values, for a total sample of 
245 sources. 



4.1.2. Lace from Lhc 

There are several empirical relations between the Lna 
and the total accretion luminosity. Since, however, we 
have an independent estimate of the Lace for some of our 
ONC sources from the 2CD, we rederive this transforma- 
tion based on our data. To this purpose, we isolate the 
sources where both Lace obtained with the 2CD and Lna 
are available. In particular, li miting to sources with at 



least a 2 cr confidence (see Sec. 3.1.3), we isolate a sam- 
ple of 148 stars. Fig. [5] shows the relation between these 
two quantities; to account for the overall uncertainties, 
including the intrinsic sca tter in the values in the Ha 
flux due to variability (e.g. Murphy et al.||20lT ), we add 
quadratically to all the estimated errors on log(Lacc/i0) 
a minimum value of 0.2, and use these combined errors 
for the linear regression. Our fit provides the following 
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Table 2 

Stellar parameters of the sourees where the 2CD method leads to a determination of Ay and h. 





RA 


DEC 




Av 


log(L./L0) 


R, 




Confidence'^ 




(J2000) 


(J2000) 


(K) 


(mag) 


(R©) 


(mag) 


fj 


20 


5:34:15.087 


-5:22:59.971 


3328 


0.56 


-1.62 


0.47 


-0.0020 


2 


30 


5:34:17.299 


-5:22:47.977 


3104 


1.09 


-1.51 


0.61 


-0.0044 


1 


47 


5:34:20.794 


-5:23:29.131 


3388 


0.51 


-1.01 


0.90 


-0.0023 


3 


50 


5:34:24.769 


-5:22:10.443 


3465 


0.22 


-0.82 


1.09 


-0.0021 


3 


63 


5:34:26.501 


-5:23:23.986 


3138 


0.95 


-0.80 


1.35 


-0.0001 


1 


69 


5:34:28.962 


-5:23:48.085 


3107 


0.42 


-1.61 


0.54 


-0.0031 


3 


70 


5:34:29.446 


-5:23:37.437 


3138 


0.08 


-1.58 


0.55 


-0.0002 


1 


77 


5:34:32.658 


-5:21:7.4565 


3632 


0.91 


-0.57 


1.31 


-0.0197 


3 


87 


5:34:30.477 


-5:21:55.949 


3124 


1.09 


-1.63 


0.52 


-0.0056 


1 


146 


5:34:48.830 


-5:23:17.984 


3705 


1.29 


-0.16 


2.03 


-0.0004 


3 


167 


5:34:50.496 


-5:23:35.331 


3435 


0.08 


-0.83 


1.09 


-0.0364 


3 



Note. — ^Orion Master Catalog entry number (Robberto et al. 2012t . *'From[Da Rio et al 



Sigma of confidence introduced m Sec. p.l.3| (This faBIe 



and references therein, 

its entirety in a machine-readable form in the online journal. A portion is shown here for 
regarding its form and content.) 



l |2012 i 

available m 
uidance 



Table 3 ^ 

Stellar parameters of the sources where we assumed Ay from |Da Rio et al.| | |2012[ | and we estimate 

Lace from the Ha luminosfty 





RA 


DEC 






log(L./L0) 


R* 




log(L^a/LQ) 




(J2000) 


(J2000) 


(K) 


(mag) 




(Ro) 


(A) 




80 


5:34:33.569 


-5:22:8.7996 


3165 


1.32 


-0.80 


1.33 


110.40 


-3.40 


101 


5:34:40.859 


-5:22:42.345 


4060 


1.15 


0.25 


2.71 


80.72 


-2.03 


104 


5:34:41.952 


-5:21:32.102 


3079 


1.24 


-1.13 


0.96 


31.52 


-4.33 


108 


5:34:41.614 


-5:23:57.485 


3106 


0.80 


-0.86 


1.29 


6.83 


-4.71 


111 


5:34:42.477 


-5:22:46.243 


3059 


1.54 


-1.13 


0.97 


70.33 


-4.00 


117 


5:34:41.816 


-5:21:49.500 


2930 


0.49 


-1.23 


0.95 


8.79 


-5.12 


130 


5:34:44.797 


-5:22:38.543 


2954 


0.83 


-1.70 


0.54 


53.61 


-4.79 


133 


5:34:46.783 


-5:21:29.147 


3495 


0.88 


-0.88 


0.99 


26.31 


-3.92 


134 


5:34:45.953 


-5:22:50.690 


3052 


2.58 


-1.11 


0.99 


11.37 


-4.78 


135 


5:34:45.987 


-5:22:47.540 


3034 


2.60 


-1.33 


0.79 


232.10 


-3.71 


144 


5:34:46.559 


-5:23:25.585 


3560 


3.40 


-0.39 


1.67 


92.06 


-2.85 



relation: 

log -r 



Note. — ''Orion Master Catalog entry number I Robberto et al.|2012 i. ^From lPa Rio et al. ( 2012 i 
and references therein. (This table is available in its entirety in a machine-readable torm m the onlme 
journal. A portion is shown here for guidance regarding its form and content.) 



(1.31 ±0.03) -log 



Herczeg & Hillenbrand 
and veiling features anc 



(2.63 ±0.13). 

(5) 

(|2008|, fitting the Balmer Jump 



correlating those estimate for 



Lace to the spectroscopical estimated E Wga, found a 
relation with similar slope (1.20 ± 0.16). |Dahm| ( |2008 1 
measured a slope of (1.18 ± 0.26) using 14 meiribers of 
the Taurus- Auriga complex, combing p revious Lace esti- 
mate w ith their spectroscopical EW Ha- |De Marchi et al. 
( 2010 1 assumed a unitary slope, implying the proportion- 
ality between Lna and L^cc, an d derived another relation 
using the same data of |Dahm] We rely on our estimate, 
and we thus use the relation in Equation [5] to obtain 
Lace for 528 additional sources with no estimate of this 
quantity from the 2CD. 

4.2. Stellar Masses and Ages 

To derive the A/ace for our sources, we still need to esti- 
mate the masses of our sources {M^, ) ; this is accomplished 
using evolutionary model interpolation of the position of 



g -2 r 




og(LH„/Le) 



Figure 5. Calibration of the Ha luminosity as an accretion indi- 
cator using the accretion luminosity (Lace) obtained with the 2CD. 
The linear regression fit shown in figure has an analytic expression 
reported in the text in Eq. [s] that permits to obtain accretion 
luminosity from the Ha luminosity estimate. 

our sources on the Hertzprung-Russel diagram (HRD). 
Moreover, this procedure gives us the age of the sources 
that will be used in the later analysis. 
In Fig. |6] we present the HRD of our sources, built 
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Table 4 

Accretion values of the sources where the 2CD method leads to a determination of Ay and La 

evolutionary models 



according to different 





log(Lacc/L0) 


(Mq) 


log Age'' 

(yr) 


log Mace'' 

(M0/yr) 


(Mq) 


log Age'' 

(yr) 


log Mace'' 

(Mo/yr) 


(Mo) 


log Age"^ 

(yr) 


log Mace"* 

(Mo/yr) 


20 


-3.26 


0.25 


7.4 


-10.40 


0.20 


7.4 


-10.30 


0.23 


7.4 


-10.36 


30 


-2.95 


0.17 


6.9 


-9.79 


0.12 


6.8 


-9.66 


0.13 


6.7 


-9.69 


47 


-2.60 


0.25 


6.5 


-9.45 


0.28 


6.8 


-9.49 


0.25 


6.5 


-9.45 


50 


-2.42 


0.29 


6.4 


-9.25 


0.33 


6.6 


-9.31 


0.31 


6.4 


-9.28 


63 


-3.59 


0.16 


6.1 


-10.07 


0.19 


6.4 


-10.15 


0.14 


5.9 


-10.00 


69 


-3.17 


0.17 


7.0 


-10.06 


0.13 


6.9 


-9.94 


0.14 


6.9 


-9.98 


70 


-4.13 


0.18 


7.0 


-11.04 


0.14 


7.0 


-10.94 


0.14 


6.8 


-10.92 


77 


-1.45 


0.40 


6.3 


-8.33 


0.42 


6.5 


-8.36 


0.43 


6.4 


-8.36 


87 


-2.99 


0.17 


7.1 


-9.91 


0.13 


7.0 


-9.77 


0.14 


6.9 


-9.83 


146 


-2.31 


0.32 


5.7 


-8.92 


0.47 


6.1 


-9.08 


0.47 


6.0 


-9.09 


167 


-1.56 


0.27 


6.4 


-8.37 


0.31 


6.6 


-8.42 


0.28 


6.4 


-8.38 



''Orion Mas ter Catalog entry number i Robberto et al.|2012|l. '"Obtained with D'Antona & Mazzitelli ( 1994 1 



Note 

models. "Obtained with Siess et al 

in its entirety in a machine-readable torni in the online journal, 
and content.) 



( 2000| models. "Obtained with |Palla fc Stahler| ( |l999| models. (This table is available 
A portion is shown here for guidance regarding its form 



Table s 

Accretion values of the sources where we assumed Ay from |Da Rio et al.| l |2012| | and we estimate Lace from the Ha 

luminosity, according to different evolutionary models 



OM'' 


log(Laec/Lo) 




log Age'' 


log Mace'' 




log Age" 


log Mace" 




log Age'' 


log Mace'' 




(Mo) 


(yr) 


(Mo/yr) 


(Mo) 


(yr) 


(Mo/yr) 


(Mo) 


(yr) 


(Mo/yr) 


80 


-1.82 


0.17 


6.1 


-8.33 


0.21 


6.4 


-8.42 


0.16 


6.0 


-8.29 


101 


-0.03 


0.40 


5.4 


-6.60 


0.75 


5.9 


-6.87 


0.80 


5.8 


-6.90 


104 


-3.04 


0.16 


6.3 


-9.66 


0.15 


6.6 


-9.63 


0.12 


6.2 


-9.54 


108 


-3.54 


0.15 


6.1 


-10.02 


0.18 


6.4 


-10.09 


0.13 


6.0 


-9.94 


111 


-2.61 


0.15 


6.3 


-9.21 


0.14 


6.6 


-9.18 


0.11 


6.2 


-9.07 


117 


-4.08 


0.12 


6.4 


-10.60 


0.10 


6.6 


-10.51 


< 0.10 






130 


-3.65 


0.11 


6.8 


-10.34 


< 0.10 






< 0.10 






133 


-2.50 


0.31 


6.5 


-9.40 


0.34 


6.7 


-9.44 


0.33 


6.6 


-9.43 


134 


-3.64 


0.15 


6.3 


-10.21 


0.14 


6.6 


-10.17 


0.11 


6.2 


-10.07 


135 


-2.22 


0.15 


6.5 


-8.91 


0.11 


6.7 


-8.76 


0.10 


6.3 


-8.74 


144 


-1.11 


0.31 


5.9 


-7.78 


0.38 


6.2 


-7.87 


0.37 


6.2 


-7.86 



Orion Ma ster Catalog entry number ( Robberto et al. l 2012||. ''Obtaine d with D'Antona & Mazzitelh 1 19941 



Note. 

models. "Obtained with Siess et al.| l |2000| | models. "Obtained with |Palla fc Stabler] | |l999| models.(This table is available 
in its entirety in a machine-readable torni in the online journal. A portion is shown here for guidance regarding its form 
and content.) 



using Tcff from Da Rio e t al. 
sity (L*) derived in Sec. [6.3 



(|2012|) and the lumino- 
a measure- 



Sources with 

ment of accretion from the 2(JD are plotted in triangles, 
whereas accretion estimates from the Ha excess are rep- 
resented by open squares. From Fig. [6] it is evident that 
our sample of Mace estimates extends well down to the 
Hydrogen-burning limit. In the very-low mass range, our 
sample includes members at all luminosities (or ages), in 
agreeme nt with the overall lumin osity spread measured 
in Orion (Da Rio et al. 2010 2012j ). We find some sources 



located on our HRU below the Main Sequence; these are 
either relatively less luminous sources seen in scattered 



light (Skemer et al. [2011 ) or non- members of the ONC, 
and will be excluded trorn our analysis. 

We assign masses and ages to our sources by inter- 
polating theoretical isochrones and evolutionary tracks_ 
on the HRD. We use m odels of D'Antona & Mazzitelli 
(pOOOl and,Palla & Stabler, ( ,1999 ) (see 

5S dist] 



(1994 
Flgl^ 



Siess et al. 
"RTF 



ed fro 



"we present the age and mass 



Tstribu- 

tions~obtained from the different evolutionary models. 
Considering only sources that appear on the HRD at 



higher than the ZAMS, our stellar sampl e includes 
730 sources w ith M, a nd age estimates from [D'Antona] 
|fc Mazzitelli] 675 using IPalla fc Stabler] models, and 697 
assuming isochrones from ] Siess et al.] 

4.2. L Completeness analysis 
The completeness of these sampl es depends on two fao 



tors: the availability of T^g from Da Rio et al.] ( ]2012 1 
and the detection of the sources in U or Ha (within the 
WFPC2 field of view). Whereas the U—h&nd photome- 
try tends to be shallow for highly reddened stars, the Ha 
photometry is much less affected by extinction, and ex- 
tends well into the brown dwarf mass range. Therefore, 
source detection is not expected to alter the represen- 
tativeness of our stellar sample. On the other hand, the 
available T^s spectral types from Da Rio et al. ( 2012 ) are 
characterized by a uniform and close to 1007o compiete- 
ness down to the H-burning limit. Therefore, we expect 
our sample of Mace to be representative of the ONC stel- 
lar population down to O.IMq, and the lower number 



of sources in our catalog compared to Da Rio et al. ( 2012 1 
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,1.50 Me D'Antona & Mazzitelli (1994) 
■ ,1 .00 M„ 



,0.10 Mj, : 
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2 - 3,50 M 



Polio & Stohler (1999) ^ 




^0 Myr - 
,190, ^1y,r: 



Figure 6. The Hcrtzsprung-Russel Diagram of our ONC sources 
with overplotted different PMS evolutionary models. Filled trian- 
gles refer to sources whose Lace has been estimated using the 2CD, 
while empty squares are sources with results derived from the Lua ■ 

is mainly due to our smaller field of view. To test this 
hypothesis, and rule out significant selection effects in 
our sample of Mace estimates, we consider full catalo g 
of available stellar parameters from 
limited to the same (smaller) field o: 



Da Rio et al. (2012), 
vIew"oFourVVTFTj2 

survey. Then we compute, as a function of mass and age, 
the fraction of these sources includ ed in our Aiacc sam- 
ple. The result, obtained assuming [D'Antona fc Mazz-| 
|itelli| models, is shown in Fig.[8j The error bars represent 
the poissonian uncertainty from the stellar numbers in 
both samples. As supposed, we do not detect significant 
trends in the fraction of accreting sources with respect 



to the stellar paramters, except for a significant lack of 
sources below 0.1 M0. Above this mass, about 2/3 of 

the ONC members have a Afacc measurement, regardless 
their mass or age. 



4.3. Mass Accretion Rates 

As show in Equation [l] Mace is derived from the mea- 
sured total accretion luminosity Lace, the stellar radius 
i?* and the mass M* . We obtain the radius of the sources 
from their Tf.g and . The obtained values of Afaec are 
reported in Table |4] and Table [5| separately for source 
with Lace estimated from the 2CD and from the Ha ex- 
cess. 

Fig. [9] shows the distribution of Mace according to 
the three evolutionary models. The distributions ap- 
pear similar to each other, spanning a range between 
~ 10"" and ~ 10-"^ Mg/yr, with a mean value of about 
r^ 1.4 • lO~^M0/yr. This is compatible with what was 
previously found for the Trapezium region ( Robberto 
et al.||2004|). 



In the same plot we also compare the distributions sep- 
arately for very low mass stars (M* < O.SM©, blue dot- 
ted histograms) and low- and intermediate mass stars 
(M* > O.3M0, red dot-dashed histograms). It is evident 
that, regardless of the evolutionary model, intermediate 
mass stars show higher mass accretion rates than lower 
mass stars. In the next section we will investigate this 
trend in more detail. 

4.4. Lace from Censored Data 
Our WFPC2 catalo g of ONC sources pr ovides pho- 



tometric upper limits (Robberto et al 



2012) for sources 



undetected in one or more photometric bands; this allows 
us to obtain, for some sources, additional upper values 
of iace- For 55 sources with available T^g and detection 
in B and /, but undetected in [/-band, we could derive 
upper limits on Lace with the 2CD method. For 8 of 
these stars with no correspondent Lho estimate, we find 
an intersection in the 2CD, leading to an upper limit on 
their Lace that can be used in our analysis. These val- 
ues are also reported in Table |4] For other 21 of these 
sources with intersection on the 2CD, Lace was already 
derived from the Lhu excess, therefore, as a sanity check, 
we compare these values with the upper limits obtained 
from U: in 14 cases (66%) we find that the upper limit 
value is higher than the derived one, while in the remain- 
ing 7 (33%) cases the two results are compatible. 

For 84 sources we have TJa-band upper limit values 
and corresponding spectral types, therefore we are able 
to derive the upper limit of Lace from Equation [5] These 
values are reported in Table [5| 

5. ANALYSIS: DEPENDENCE OF Lacc AND Mace ON 
STELLAR PARAMETERS 



5.1. L, 



ace L>^ 



In Fig. 10 we plot Lacc/L© as a function of L^/Lq. We 



stress that these two quantities are independent of the as 
sumed evolutionary model, allowing a less biased analysis 
of the dependence of accretion with the stellar properties. 
Because of this, such rela tion has been long investi gated: 
Clarke fc Pringle ( 2006 1 , reporting results from 'Natta 
et al.~( |20U6 ), showed an almost uniform distribution of 
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1.2 
1 .0 
0.8 
0.6 
0.4 

0.2 
0.0 



-0.5 
log M (Mg) 



0.0 



Figure 8. Frac tion of the sources in our samples witii respect 
to the sample of |Da Rio et alTj (|2012[ l, that is complete at 100% 
down to the H-burning hmit. We see that we have a ^70% of 
completeness down to the H-burning limit. 

sources on this plot between the upper bound Lace ~ 
and the lower observational detecti on threshold, with a 



relation Lace oc L% with {3 = 1.6. [Tilling et al | ( 2008 1 
were able to theoretically reproduce the observed disper- 
sion for > 0.4 Mq, even if they could not reach the 
lowest values of Lace de tected observationally. Assuming 
D'Antona & Mazzitelli evolutionary mo dels, they used a 
depen dence of Mace cx t~^' with rj = 1.5 (Hartmann et al 
1998[) obtaining ^ = 1.7. 



Pig. 10 shows that our plot is widely populated between 
Lace L* and Lace 0.01 L*, with several sources below 
this rang e, down to Lace ~ 0.001 L*. The upper locus 
found by Natta et al. (2006) is still present at Lacr ^ L* 



in our sample. We observe 21 sources over this level, but 
17 of which are compatible within 2 sigma to L^^c ^ L^ . 



The remaining 4 sourc es may be " continuum stars" ( Cal- 
vet & Gullbring 1998), i.e. sources where the accretion 
component dominates over the stellar emission. If this 
were the case, their spectral types might also be highly 
uncertain, with also an impossibility to determine the 
correct L*. 

We note also that the density of sources with L* > 
Lq is much lower than that of lower luminosity sources, 
probably because of observational selection effects. In 
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Figure 9. Distribution of mass accretion rates according to the 
three different evolutionary models. In each plot we also present 
the partial distributions separating the sample in two mass bins, 
highlighting the different averages. 

particular, as will be pointed out in Sec. |5.2.3[ for higher 
mass stars our 2CD does not allow us to detect low values 
of Lace- Viceversa, the lack of high values of Lace for 
L.^, < 0.3 Lq occurs in a region widely populated by 
upper limit detections, and thus cannot be attributed to 
observational limits. 
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Log(L,/Le) 

Figure 10. Accretion luminosity vs stellar luminosity. Values 
obtained from the 2CD are plotted as triangles, while squares rep- 
resent data obtained from Ha excess. Dot - dashed lines represent 
different ratio between accretion and stellar luminosity and the 
solid line is our best fit (slope~1.68) from a linear least-squares 
regression which accounts for errors on both L» and Lace- The fit 
assuming only results from the 2CD is showed with a dotted line, 
while the one from only Ha luminosities with a dashed line. 

Fitting our results with a hnear least-squares proce- 
dure, accounting for the errors (corrected as in Sec. 3.1.3 
for both stellar and accretion luminosity), we obt ain a 
slope f3 ~ 1.68 ± 0.02. This is also shown in Figure 10 



as 



a solid line. This v alue is compa t ible w ith those obtained 
obser vationally by Natta et al. (2006) and theoretically 

that the slope 



from 'Tilling et al. 



We note. 



slightly changes w. 
determined from the 2CD 
(dotted line), or from Lna, 
(dashed line). 



however, 

len selecting only sources with Lace 
for which /3 = 1.73 ± 0.02 
for which /3 = 1.59 ± 0.04 



5.2. Mace vs Age and 
5.2.1. Fit to the data 

To investigate the variations of Aface with respect to 
the stellar age and mass, we consider all these 3 quanti- 
ties in a single a 3D space. To this purpose, we neglect 
sources with suspicious ages, specifically, too young (log 
(t/yr) < 5.5) or too old (log (t/yr) > 7.3). 

We perform a minimum Xicd accounting for errors 
on the single values, assuming two different forms: a) a 
simple plane (logMacc= A + B ■ log t + C ■ log A/* ) and b) 
a warped surface (logAfaec= A + B ■ log t + C ■ log Af* -f 
D ■ log t ■ log ) . Specifically, the first form is the sim- 
plest function to describe correlations of Aface with the 
stellar parameters, the second allows for a mixed term 
that accounts for a possible different evolution timescale 
of Aface with stellar mass. To perform the fit we use a 
standard Levenberg-Marquard non-linear regression. In 
the two cases the resulting x^cd ^^'^ comparable (^ 7—9), 
and the second one being slightly smaller. This decrease 
of Xrcd is '^ot surprising, given that the first function is 
a subcase of the second, the latter having an additional 
parameter. Therefore, in order to quantitatively decide 
whether the most complex model is truly more repre- 
sentative of the real distribution of the data, or whether 
the simplest model is good enough given the uncertain- 



.4 r 
^6 ' 




11 and described by the 



Figure 11. Snapshot of our 3D diagram of Mace with respect to 
both M, and age, together with our best fit surface (see text). The 
analytical expression of the latter is also reported in Eq.|6]and the 
coefficients of the fit are reported in T able[6| The full 3D animation 
is available electronically. In Fig. [12] different views of this surface 
are shown. 

ties, we run a statistical F test. This test compares the 
relative increase in the residual sum of squares obtained 
by reducing the complexity of the model, with the rela- 
tive increase of the degrees of freedom. Then, the test 
derives the probability that the simpler model is suit- 
able enough. In any case, we find this probability lower 
than 0.1%, meaning that introducing the fourth param- 
eter D in our model (i.e., a warped surface instead of 
a flat plane) leads to a more representative match with 
our data, and this is not a statistical casualness. Our 
fitted surface is shown in Fig. 
following equation: 



log Afaec = A + B ■ \ogt + C ■ logAf* + D ■ logi • log A<^) 

where the parameters' values are reported in Table [6] also 
in the cases where we consider only sources with results 
obtained from the 2CD or with accretion estimated using 
the Ha luminosity. 

These relations imply that both parameters play to- 
gether a role in the evolution of mass accretion. To bet- 
ter understand this behavior, and in particular the role 
of the "mixed" term described by the parameter D, we 
separately investigate the dependence of Aface on age or 
mass; also, in every case, we divide the sample in diffe- 
rent ranges for these parameters. This is illustrated in 
Table[7]for the dependence on age, and in Table[8]for the 
dependence on mass, and the plots are shown in Fig. |12[ 

5.2.2. Discussion 

If we consider only the dependence of Aface on stellar 
ages (AfaccOt t~^), om fit implies that a si milar slope as 
that proposed by Hartmann et al. ( 1998| ) {rj — 1.5) is 
obtained for Af* ~ 0.13 Mq. b'or higher mass stars we 
find decreasingly lower values of rj with mass. Table [7] 
in particular, shows that the values of rj decrease mono- 
tonically regardless of considering all sources, only those 
with results obtained from the 2CD or with the Ha lu- 
minosity. Thus, our data indicate that there are different 
accretion evolutionary timescales for sources with diffe- 
rent mass. Surprisingly, for sources with Af* > 0.5AfQ 
the slope is smaller than 1. This value is not compatible 
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Table 6 

Coefficients of Eq. IgI obtained fitting the values of Mace as a function of both Af, and 
the age of the sources for different subsamples 



Sample 


A 


B 


C 


D 


All sources 
Only 2CD results 
Only Ha results 


-(5.12 ±0.86) 
-(4.70 ± 1.29) 
-(3.29 ± 1.36) 


-(0.46 ±0.13) 
-(0.56 ±0.21) 
-(0.71 ± 0.21) 


-(5.75 ± 1.47) 
(0.35 ± 0.35) 
(0.93 ± 0.34) 


(1.17 ±0.23) 
-(0.54 ± 2.20) 
-(4.01 ± 2.19) 



Note. — The sample is divided in subsamples according to the method used to obtain 
Lace, and thus A/acc. 
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with the similarity solution model framework (Lynden- 
Bell &: Pringle||1974[ |Hartmann et al.||1998[ ). Indeed, 
it we suppose that the disk evolves with time follow- 
ing a simple power law (Afdisi< oc t^"), we have that 
M = dMd/dt (X thus 77 = a + 1 > 1 in order to 

have a disk mass decreasing with time. 

On the other hand, according to the self-similar mod- 
els. Mace scales with time as a simple power law only 
in the asymptotic regime, when the age of the system 
t ^ t^, where ti, is the viscous timescale of the disk 
(set by the initial conditions and by the relevant vis- 
cosity). At earlier times, the relationship between Mace 



Figure 12. Mace evolution as a function of individual age and mass of the sources for different ranges of these parameters. The red line 
is plotted according to the value of the fit obtained using all the sources in the sample (Eq. [6| and for the mean value of the parameter in 
the range written on the picture. 

and t is indeed shallower. Our results might thus be ex- 
plained within the framework of the self-similar solutions 
if higher mass stars are not yet in the asymptotic regime, 
implying that the viscous time is a growing factor of the 
stellar mass. Note that this trend is t he opposite of what 
suggested by Alexander & Armitage ( [200 6 ) in the brown 



dwarf regime to explain the M*-Macc correlation. 

Alternatively, one might question the validity of the 
self-similar models themselves in describing, from a sta- 
tistical point of view, the evolution of protostellar disks, 
and in particular the assumption that viscosity scales 
with radius as a simple power law. Indeed, different 
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Table 7 

Afacc VS. ago relation coefficients 



ri(2CD) 



O.13M0 

0.2Mq 

0.3Mq 

O.5M0 

O.8M0 

IA/q 



1.50 ± 0.26 
1.28 ±0.21 
1.07 ± 0.18 
0.81 ± 0.15 
0.58 ±0.13 
0.46 ±0.13 



0.87 ± 0.41 
0.80 ± 0.32 
0.74 ± 0.28 
0.66 ±0.23 
0.59 ±0.21 
0.56 ±0.21 



1.53 ±0.40 
1.36 ±0.32 
1.20 ±0.28 
0.99 ±0.23 
0.80 ±0.21 
0.71 ±0.21 



Note. — Results of the fit A/acc obtained consid- 

ering all the sources (second column), sources with results 
obtained with the 2CD (third column) and with the Ha lu- 
minosity (fourth column). Values of the slope are different 
for stars belonging to subsamples with different mean 
(first column). 



Table 8 

Afacc VS. Mt relation coefficients 







b 


fe(2CD) 


b{Ha) 


age r- 


- 0.8 Myr 


1.15 ± 2.00 


1.51 ± 3.02 


1.50 ± 2.97 


age r- 


- 1 Myr 


1.26 ± 2.02 


1.54 ± 3.04 


1.59 ± 2.99 


age 


- 2 Myr 


1.61 ± 2.06 


1.65 ±3.12 


1.87 ± 3.06 


age r- 


- 5 Myr 


2.08 ± 2.13 


1.79 ± 3.22 


2.24 ± 3.16 


age r- 


- 8 Myr 


2.32 ± 2.16 


1.86 ±3.27 


2.43 ± 3.21 


age 


- 10 Myr 


2.43 ± 2.18 


1.89 ±3.29 


2.52 ± 3.23 



Note. — Results of the fit Mace obtained consid- 

ering all the sources (second column), sources with results 
obtained with the 2CD (third column) and with the Ha lu- 
minosity (fourth column). Values of the slope are different 
for stars belonging to subsamples with different mean age 
(first column). 

physical processes are expected to redistribute angular 
momentum within the disk at different radii, with the 



magne to- rotational instability (MRI, Balbus & Hawley 



(1992)) dominating at s mall distances from th e star, the 



gravitational instability (Lodato & Rice 



20041 dominat- 



ing in the colder outer parts, and with the p ossibility of 
having extended dead zones ( Gammie||1996 ) in between 
the two. 

For what concerns the dependence of Mace over mass 
(see Table [s]) modeled as MaccOc A'/^, the trend is similar 

regardless on the method used to estimate Mace- In par- 
ticular, the slope in higher, thus the relation is steeper, 
when considering older ages. This also shows that at 
older ages more massive sources maintain an higher ac- 
cretio n ra te. 

12 illustrates once more these trends. We plot the 
values against age and for diffe- 
In this figure the slopes of 



Fig._ 
evolution of our M. 
rent ranges of the parameters 

the overplotted lines are obtained from Equation |6j sub 
stituting to the projected parameter (mass for the left 
panels, age for the right ones), the mean value within 
the subsample used for each particular panel. 

The slower temporal decay of Mace we measure in the 
intermediate-mass range {M^ > 0.4 Mq) is in agreement 
with estimates of mass accretion rates fr om Ha photom- 
etry o btain ed in the Magellani c Clouds ( jPe Marchi et al. 
' 2010[ 2011 Spezzi et aL|2012| and in the massive galac- 



tic cluster In G(J3fci03 (Beccari et ai. 2010). Indeed, in 



all these studies, the investigated stellar samples do not 
reach the low-mass regime, due to poor photometric sen- 
sitivity. This suggest, as inferred by the authors of these 



works, that their relatively high values of Mace measured 
for isochronal ages of '-^lOMyr compared to galactic low- 
mass stars indicate a selection effect because of higher 
masses. 

5.2.3. Caveats 

Our derivation of A^acc using the 2CD may be, in prin- 
ciple, affected by some modest selection effects. This, 
as mentioned, is located in the fact that for high Tes 
sources, the color displacement due to accretion is very 
small, limiting the possibility to identify weak accretors 
among intermediate mass stars. Specifically, Fig.[3]shows 
that for [B — I) < 2, the 1% iacc/-^tot line is undistin- 
guishable from the zero accretion one. This is simply 
because stars with high Tofi emit a non-negligible frac- 
tion of their flux in the C/-band, "hiding" a possible small 
accretion excess. As mentioned, however, this problem 
is highly mitigated since we have also derived Lace from 
the Lna excess for all the sources for which Mace is not 
derived from the 2CD. In the end, with the Ha mea- 
surements, and even considering upper limits (see, e.g. 
Fig. 10) the relations we find do not appear to suffer 
from this mentioned selection effect. 

As already introduced in Sec. |4.2[ the accuracy of the 
isochronal ages to re present the true stella r age s has been 
recently questio ned. Baraffe et al. (2009) and Baraffe & 
Chabrier ( 2010[) claimed that the Doserved spread of 
m the MiId of star forming regions can be explained 
as a consequence of the protostellar phase accretion his- 
tory. In particular, episodic events of intense accretion 
(Macc^^ lO~'*M0/yr) during this phase can produce a 
large spread of values of stellar radii for stars of identical 
age and mass. Their results are based on the assump- 
tion of "cold" accretion and initial masses Mj between 



1 M,7 and 0.1 Mq. On the other hand, Hosokawa et al. 
(2010) demonstrate that the apparent spread produced 
by different accretion histories may be important only 
for sources with T^ff ^ 3500 K (^ 30% of our sources). 
In their analysis they assume a fixed — 0.01 Mq, and 
this could suggest that the results obtained should be re- 
vised changing the value of this parameter. Thus, from a 
theoretical point of view the debate is still ongoing and 
a unified picture is not yet reached. At the present mo- 
ment, therefore, we can only use the age inf erred from 



the evolution a.ry tracks used in the past (e.g. D'Antona 
& MazzitelliJ ([l994)), waiting for new isochrones based 



on episodic accretion models. 

We stress that if the isochronal ages we assume are 
inaccurate, because of the aforementioned effects, this 
results in an overall "horizontal" broadening of our re- 
sults in the Mace vs. age plane, but the overall trends we 
measure cannot be strongly influcenced. 



5.3. 



Mace vs position 



We look for spatial variations of Afacc within the ONC. 
Such variations, in fact, might imply that the accretion 
is subject to environmental effects, such as the infiuence 
of the massive stars present in the central region of the 
cluster, in particular 9^ Orionis C, an 0-type star. 

In order to do that, we divide our sample in 5 subsam- 
ples, each with the same amount of sources, according 
to the distance of each object from 9^ C. We show with 
black squared symbols in Fig. 13 the mean values of the 
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B'C (p= 



Figure 13. Mean values of Mace at increasing distances from 
6^ Orionis C. The black points represent tlie values using all the 
sources in the sample and dividing them in 5 subsamples of equal 
number of sources. Red points represent the mean values when se- 
lecting for each subsample only the sources with M» > 0.3Mq and 
5.5< log(age/yr)<7.3. In the first case a trend of lower Mace at 
higher distances is clear, but can be due, at least partially, to selec- 
tion effects, as we show with the red points, where the distribution 
is flatter. 

-M^acc plotted against the mean value of the distance for 
each subsample. It is evident that at smaller distances 
from the central star sources appear to have an higher 
mean Mace that decreases monotonically with the dis- 
tance. This result, however, does not necessarily indi- 
cate a real trend. In fact, around the Trapezium cluster, 
the detection limit is generally poorer than at larger dis- 
tances, due to the higher nebular background. Therefore, 
low-mass stars (whose Mace are generally lower, as de- 
scribed in the previous sections) are less easily detected. 
Therefore, we select in each subsample only the sources 
with < O.3M0 and with 5.5 < log(t/yr) < 7.3, to 
have similar selection effect in e ach bin. The result is 
shown with the red circles in Fig. [l3| the values of Mace 
as a function of distance have a much flatter distribu- 
tion. Only at distances > 1.5 pc the mean value of Mace 
is lower. 

6. CONCLUSION 

We have presented a study of mass accretion rates 
(Mace) in the Orion Nebula Cluster, for an unprecedent- 
edly large (~ 700 stars) sample of PMS stars. This al- 
lowed us to perform a thorough statistical analysis of 
the dependence of this quantity on the central stellar pa- 
rameters, and investigate this phase of the stellar mass 
build-up and disc evolution. 

We based our study on IIST/WFPC2 photometric 

data over a large field of view, and derived Mace using 
two different accretion tracers: the {U — B) excess and 
the Lhu- We study the systematic dependence of Aface 
with respect the age of the sources, their mass and posi- 
tion within the region. In particular, the Mace is found 
to vary with age and mass altogether. Our final relation 
between all these quantities is given by Equation|6]and in 
Table § and our best fit relation using all the sources in 
the sample is given by: log(Macc/M0-yr)=(-5.12 ± 0.86) 
-(0.46 ± 0.13) •log(i/yr) -(5.75 ± 1.47)- log(M*/M0) + 
(1.17 ± 0.23)-log(</yr) • log(M*/M0). We clearly find 
that the Mace increases with stellar mass, and decreases 



over evolutionary time. 

Interestingly, we also find evidence that for more mas- 
sive stars the decay of Mace with time is much slower than 
for lower stellar masses. Similarly, for older stars, the de- 
pendence of Mace with appears significantly steeper. 
This might imply that these objects are not in the asymp- 
totic regime (i.e. when t t^, where Mace oc t^**), or 
that the hypothesis of a simple dependence of the viscos- 
ity v oc is probably not compatible with our observa- 
tions. In particular, we suggest significant discrepancies 
of our results with respect to the self-similar parametriza- 
tion for sources with masses higher than ~ O.5M0. 

This work was made possible by GO program 10246 
of the Hubble Space Telescope, which is operated by the 
Space Telescope Science Institute. 

C.F.M. would like to thank James Pringle and Philip 
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Figure 14. Bandpass of the filter F336W/WFPC2. Note tlie red leak component at ^ 0.75^m, where the filter has a window of 
transmission with 1% of the peak transmission at ~ 0.35^m. 



APPENDIX 



RED LEAK 

As mentioned in Sec. |2.1[ the filter F336W of WFPC2 is affected by red leak. This is evident from Fig. [T4j where 
the band profile is shown. At A ~7300 A, the filter shows a transmission window whose peak throughput is ~0.5% 
that in the [/-band wavelength range. There is also a third peak at '-^ 1 /im, but this is completely negligible, since its 
bandpass is less than 10^^ of the main peak and the detector Quantum Efficiency is nearly zero. 

The presence of red leak in F336W does not significantly affect the stellar photometry of blue sources, as most of 
the flux is emitted at short wavelengths. For red objects, however, the leak contribution may be significant and, in 
extreme cases, even dominant. In order to accurately account for the J7-band fluxes of our ONC sources, and prevent 
introducing a systematic overestimate of the accretion luminosity (iacc) and mass accretion rate (Mace), we must 
estimate and correct for the spurious red-leak excess. 

The number of photons collected in the leaking part of the filter depends solely on the observed flux at the wavelength 
of the leak (referred to as J/icak, in magnitudes). To estimate this flux, we considered our WFPC2 J-band photometry 
and computed the color term (i/ieak — I) by means of synthetic photometry. As this col or term depend s mainly on Tcs 



(2012), and assuming 



and Ay, we computed it for each source conside ring the stellar parameters from Da Rio et al. 
the atmosphere models from AUard et al. (2010), corrected as explained in Appendix [b| 

We then converted t/ioak from magnitudes to photons/second. To this end, one must derive the zero-point of the 
leak window of the filter, i.e, the number of photons passing through the lea k for a magnitude m=0. This zero point 



is simply the zero-point of the whole F336W filter {ZP p^^Qw , |Dolphin||2009 ), scaled to a term equal to the fra ction of 
U -ban d photons that pass through the leak. We derived this latter term from the calibrated spectrum of Vega of |Bohliii| 
(2007) and the measured F336W throughput. We find that the fraction of photons in the red leak part of the spectrum 
of Vega is 0.551% of the total number of [/—band photons. Therefore, we obtain countsicak = lO""-^'^'^'"''''^'^^^ -0.00551. 

Finally, we considered our observed U magnitudes, converted them in units of counts (using again "ZPf^^qw), 
subtracted the leak contribution, and transformed them once again in magnitudes. The result is the leak-corrected 
[/—band magnitude, or C/noicak- The leak offset in magnitude, for our stellar sample, has a mean value of about 0.5 
mag. We refer in all the text to the U band as the corrected one, using the extended notation C/noicak only when 
needed, for clarity. 

ISOCHRONE CALIBRATION 

In order to calibrat e the approp r iate is ochrone for our sample, we start by considering a family of synthetic spectra, 
the BT-SETTL from AUard et al. (2010). As shown in their work, these models seem to reproduce the observed colors 
of stars and brown dwarfs in the J and K band. However, their validity in the optical range has not been dem onstrated. 
We therefore test their accuracy taking advantage of the J band photometry on the ONC from Robberto et al. (2010). 
The concept is the following: assuming the validity of the predicted J-band colors as a function of I'es, it the predicted 
/-band colors is also correct, then the synthetic (/ — J) should also agree with the data. Otherwise, a correction to 
the synthetic /-band magnitude is required. Subsequently, the same procedure can iterated to the B and U bands as 
well. 

Fig. 15 illustrates this test: the upper panel shows the obser ved (grey open circ les) and extinction-corrected (black 



accretors (Da Rio et al. 



these assumptions, the black 



2010 



dots, 



dots) pEot ometry of our ONC memb ers, assuming the Ay from|Da Rio et al 

2009), i.e. sources with log(L; 
therefore 



( |2012[ ) , and limited to the candidate weak 
1.5 and EW^Q^Daffioog < 10 A. Under 



/itot) DaRiolO < 

trace the photospheric colors of O NC sources of differen t temperatures 



The red dashed line represents a 1 Myr synthetic isochrone from the models of Baraffe et al. (19981, converted into 
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Figure 15. Procedure used to calibrate the models by AUard et al.|(|2010[ l in our photometric bands. The red dashed line represents the 
color computed trough synthetic photometry o n the BT-aJU 1 "1 L mo dels. Urey open circles are the observed colors of the sources, grey filled 
one the corrected values according to Ay from |Da Rio et al.| ( |2012| and black filled points are corrected values of low accretors according 
to |Da Rio et aL| ( |2010'l l. We use the latter ones to determine tne ernpirical calibration, that is the blue line. 

colors assuming the BT-SETTL atmosphere models models. We find that the predicted colors are systematically 
underestimated for Toff < 3200 K. We thus interpolate the empirical I — J color locus (blue solid line), and determine 
the offset between the two. We calibrated the /-band magnitudes of the isochrone as a function of Toff, by applying 
this correction to the synthetic J— band magnitudes. 

From the calibrat ed / -band photospheric magnitudes, we iterate the same procedure first to the B — I colors, in order 
to calibrate B (Fig. 15 middle panel), then to the U — B color, to calibrate the intrinsic U magnitude (lower panel). 
The selection of the candidate non-accretors is particulary critical for these two bands, which are highly influenced by 
accretion excesses. 

We stress that for practical reasons we performed a calibration of the synthetic isochrone in terms of magnitudes, but 
for our purposes (the analysis of our 2CD) only the accuracy of the colors as a function of temperature are relevant. 
Also, colors are independent on distance or small differences in stellar radii, therefore a luminosity spread for ONC 
members, due to e.g., a spread in stellar ages, does not affect our method. We recap the empirical estimate of the 
WFPC2 colors found for the cahbration of the BT-SETTL models in Tabled 



